Power-law distributions in economics: a nonextensive statistical approach 
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The cornerstone of Boltzmann-Gibbs (BG) statistical mechanics is the Boltzmann-Gibbs-Jaynes- 
Shannon entropy Sbg = ~k J dx f(x) In f(x), where k is a positive constant and f(x) a probability 
density function. This theory has exibited, along more than one century, great success in the treat- 
ment of systems where short spatio/temporal correlations dominate. There are, however, anomalous 
natural and artificial systems that violate the basic requirements for its applicability. Different phys- 
ical entropies, other than the standard one, appear to be necessary in order to satisfactorily deal 
with such anomalies. One of such entropies is S q = k (1 — j dx [f(x)] q )/(l — q) (with 5*i = Sbg), 
where the entropic index q is a real parameter. It has been proposed as the basis for a general- 
ization, referred to as nonextensive statistical mechanics, of the BG theory. S q shares with Sbg 
four remarkable properties, namely concavity (Vq > 0), Lesche-stability (Vq > 0), fimteness of the 
entropy production per unit time (q € 5R), and additivity (for at least a compact support of q in- 
cluding q = 1). The simultaneous validity of these properties suggests that S q is appropriate for 
bridging, at a macroscopic level, with classical thermodynamics itself. In the same natural way 
that exponential probability functions arise in the standard context, power- law tailed distributions, 
even with exponents out of the Levy range, arise in the nonextensive framework. In this review, we 
intend to show that many processes of interest in economy, for which fat-tailed probability functions 
are empirically observed, can be described in terms of the statistical mechanisms that underly the 
nonextensive theory. 
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I. INTRODUCTION 

The concept of "entropy" (from the greek Tpown, transformation), was introduced in 1865 by Rudolf Julius Em- 
manuel Clausius in the context of Thermodynamics [lj. This was motivated by his studies on reversible and irreversible 
transformations, as a measure of the amount of energy in a physical system, that cannot be used to perform work. 
More specifically, Clausius defined change in entropy of a thermodynamic system, during some reversible process 
where a certain amount of heat SQ is transported at constant temperature T, as 5S = 5Q/T. We can consider en- 
tropy as the cornerstone of Thermodynamics, since all the thermodynamical principles involve, directly or indirectly, 
this fundamental concept. 

The first connection between the macroscopic Clausius' entropy of a system and its microscopic configurations was 
done by Ludwig Boltzmann in 1877 Studying the approach to equilibrium of an "ideal" gas Q, he realized that 
the entropy could be related to the number of possible microstates compatible with the thermodynamic properties of 
the gas. For an isolated system in its terminal stationary state (thermal equilibrium), Boltzmann observation can be 
expressed as 

S BG = k\nW, (1) 

where A: is a positive constant and W the number of microstates consistent with the macroscopic state. This famous 
equation, known as Boltzmann principle, can be regarded as a pillar in the foundations of statistical mechanics, which 
aims to describe the thermodynamic state of a system through the statistical approach of its microstates. 
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For systems which are not isolated but instead are in contact with some reservoir (of heat, particles, etc.), it is 
possible to derive (under some assumptions), from Eq. £Q) Q , the celebrated Boltzmann-Gibbs entropy 

w 

Sbg = -k^Pi In pj, (2) 

i=l 

where pi, subject to the normalization condition Yl^L-^ Pj = 1, is the probability of the microscopic configuration 
i. In particular, for equiprobability, pi — 1/W, V^, then Boltzmann-Gibbs entropy J5J) reduces to QJ. Since it 
refers to microscopic states, the Boltzmann principle should be derivable from microscopic dynamics. However, the 
implementation of such calculation has not been yet achieved. Consequently, BG statistical mechanics still remains 
based on hypothesis such as Boltzmann's Stosszahlansatz (molecular chaos hypothesis) and ergodicity 0. It can 
be easily shown that entropy ijSJ) is nonnegative, concave, experimentally robust (or Lesche- stable) pj,Hj|, and leads to a 
finite entropy production per unit time 6]. Moreover, it is additive. In other words, if A and B are two probabilistically 
independent subsystems, i.e., pfj +B — pfpf ' •> then it is straightforwardly verified that 

Sbg {A + B) = S bg (A) + S BG (B) , 

hence, if we have N subsystems, Sbg{N) = NSbg(^), where the notation is self-explanatory. More generally, when 
correlations are "weak" enough, Sbg is extensive, i.e., such that the limjv^oo Sbg(N)/N is finite. 

Despite the lack of first-principle derivations, Boltzmann-Gibbs statistical mechanics has a history plenty of successes 
in the treatment of systems where short spatio/temporal interactions dominate. This kind of interactions favor 
ergodicity and independence properties, necessary in Khinchin's approach of Sbg E9- Consequently, it is perfectly 
plausible that physical entropies other than the Boltzmann-Gibbs one can be defined in order to suitably treat 
anomalous systems, for which ergodicity and/or independence are not verified. As examples of anomalies we can 
mention: metastable states in long-range interacting Hamiltonian systems, metaequilibrium states in small systems 
(i.e., systems whose number of particles is much smaller than Avogrado's number), glassy systems, some types of 
dissipative systems, systems that in some way violate ergodicity, and, last but not least, systems with non-Markovian 
memory, like it seems to be the case of financial ones. Generically speaking, systems that may possibly have a 
multifractal, scale-free or hierarchical structure in the occupancy of their phase space. 

Motivated by this scenario, one of us proposed in 1988 the entropyQ 

w 

S q = k '- 1 ( ? 6»), (3) 

i - q_ 

which generalizes Sbg, such that lim g _,i S q — Sbg, as the basis of a possible generalization of Boltzmann-Gibbs 
statistical mechanics 8] and where the entropic index q should be determined a priori from microscopic dynamics. 
Just like Sbg, S q is nonnegative, concave, experimentally robust (or Lesche- stable) (Vg > 0), and leads to a finite 
entropy production per unit time (|9j and references therein). Moreover, it has been recently shown 10] that it is also 
additive, i.e., 

N 

S q (A 1 + A 2 + ... + A N ) = J2 S q (Ai) , 

i=l 

for special kinds of correlated systems, more precisley when phase-space is occupied in a scale-invariant form. 

Since its proposal, entropy © has been the source of numerous results in both fundamental and applied physics as 
well as in other scientific areas such as biology, chemistry, economics, geophysics and medicine [ill ]. It is our purpose 
here to review some new results concerning applications to economics and more specifically to the theory of finances. 



II. PROBABILITY DENSITY FUNCTIONS FROM THE VARIATIONAL PRINCIPLE 



Before discussing some specific quantities of financial interest, let us see the form of the probability density function 
naturally derived from the variational principle related to entropy |3"|. We consider its continuous version, i.e., 
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corresponding to the generalized mean and variance of a relevant quantity x, respectively. It is noteworthy that, 
averages weighted with a function of the probabilities as in (0-© allow to mimic the way individuals behave in face 
to risky choices^^ (see also [T^|). In fact, the prospect theory proposed by Kahneman and Tverskv|l4j for analyzing 
decision-making under risk, is founded on the concept of "decision weights" that can be modeled by akin functional 
forms 15j. 

From the variational problem for J3J under the above constraints, one obtains 



p(x)=Ag l + (q — l)Bg(x — fl q y 



(q < 3) ; 
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where, 
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and 



B q = [(3 - 1) a*] . 

The upper bound q — 3 guarantees normalizability. Defining the q- exponential function as 



[1 + (1 - g) x} 1 -" 



{el 



{e x q = if 1 + (1 — q)x < 0) we can rewrite the probability density J7J) as 



p (x) = A q 



(8) 



hereon referred to as q- Gaussian probability density function. 

For q = §4^, the q-Gaussian form recovers the Student's i-distribution with m degrees of freedom (m = 1,2,3,...). 
Consistently, for q > 1, JHJ presents an asymptotic power-law behavior. Also, if q = with n = 3,4, 5, . . ., p{x) 
recovers the r-distribution with n degrees of freedom. Consistently, for q < 1, p(x) has a compact support which is 



defined by the condition \x — fj, q \ < 

Many other entropic forms have been introduced in the literature for various interesting purposes. One of the 
advantages of entropy © is that it yields power-law tails, which play a particularly relevant role, as well known. 

Let us recall succinctly the two basic central limit theorems: (1) A convoluted distribution with a finite second 
moment approaches, in the limit of N — > oo convolutions, a Gaussian attractor; (2) A convoluted distribution with a 
divergent second moment, approaches, in the same limit, a Levy distribution L 1 (x) (with < 7 < 2) 16]. However, 
through dynamics different from the convolution one, for instance with some sort of memory across successive step s 
(i.e., nonindependence of the successive distributions), different ubiquitous distributions might arise (see also |10|). 
If the Laplacian term in the linear diffusion equation is a standard second derivative or a fractional derivative, the 
Gaussian or the Levy distributions are respectively attained as solutions. However, much more complex and rich 
dynamics clearly exist in nature, for example those associated with a variety of nonlinear Fokker-Planck equations 
involving nontrivial correlations, multiplicative noise or other effects, as we will see in the following Section. Moreover, 
simple convolutions allow only for asymptotic behavior like the q = 1 (Gaussian) or q > 5/3 (Levy distributions) 
ones. But they do not allow fat-tailed distributions associated with 1 < q < 5/3. However, many complex systems in 
nature as well as in social sciences exhibit exponents which precisely belong to that interval. An example of financial 
interest is exhibited in Fig^ 
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FIG. 1: Probability density function of returns, P(r), versus return, r. The symbols represent P(r) for the Dow Jones Industrial 
daily return index from 1900 until 2003. The solid line represents the best q-Gaussian numerical adjust with q = 1.54 and 
a\ = 0.338 (as obtained in [l7|) and the dashed line a Gaussian fit. 



III. UNDERLYING STOCHASTIC PROCESSES 



The Gaussian distribution, recovered in the limit q — ► 1 of expression JSJ), can be derived on a variety of grounds. 
For instance, it has been derived, through arguments based on the dynamics, by L. Bachelier in his 1900 work on price 
changes in Paris stock market, and also by A. Einstein in his 1905 article on Brownian motion. In particular, starting 
from a Langevin dynamics, one is able to write the corresponding Fokker-Planck equation and, from it, to obtain 
as solution the Gaussian distribution. Analogously, it is also possible, from certain classes of stochastic differential 
equations and their associated Fokker-Planck equations, to obtain the distribution given by (JBJ- 

In this section we will discuss dynamical mechanisms which lead to probability functions with asymptotic power-law 
behavior of the g-Gaussian form. 



A. Stochastic processes with additive multiplicative structure 

Microscopic dynamics containing multiplicative noise may be encountered in many dynamical processes and, due to 
its significance, has been the subject of numerous studies in the last decades. The presence of additive noise, besides 
the multiplicative contribution, is in fact a quite realistic ingredient since not all the fluctuations are processed 
multiplicatively. In previous work0, we considered a class of differential stochastic equations of the form 



x = f(x)+g(x)((t) + n(t), (9) 

where /, g are arbitrary functions of the stochastic variable x, and (^(t),r](t), two independent zero-mean Gaussian 
white noises with variance M 2 and A 2 , respectively. We have shown that for special forms of the deterministic force, 
namely f(x) — — r )g(x)g' (x), the stationary probability density functions (by using the ltd prescription) are of the 
form 



P. (a) oc e -«*)] 2 , (10) 



where f3 = (7 + M 2 )/A 2 and 



1 + 2Af 2 / 7 
1 + M 2 h 



(11) 



From the point of view of entropy S q , the density function 110(1 derives from the variational principle, under the 
constraints of normalization and generalized variance of g(u). 

In particular, g-Gaussian distributions can be derived from a stochastic process of the linear form[ll [3 HJ UK 



x = -jx + x((t) + r](t) . 



(12) 



In fact, its associated Fokker-Planck equation is 



DP _ d(xP) A 2 d 2 [{l + (M/A) 2 x 2 )P] 

dt ~ 1 dx + 2 dx 2 ( ' 



that has the alternative form 



dP _ d{xP) D dP u 

~dt _7_ foT + 2" "fe 5 "' (1 ' 

where v = 2 — q, and D is a constant related to the other model parameters. Eq. I|14|l is a nonlinear diffusion equation, 
familiarly known as porous media equation. Its steady state solution has the form 



P s (x) oc e~^ 2 (15) 

with q and (3 defined as above. In the particular case v = 1 = q the standard Gaussian steady state is obtained, 
corresponding to a purely additive process. 

Taking into account that empirical returns where found to follow a (/-Gaussian distribution (see also Fig. 
Eq. I|12(l . complemented by the Ito prescription, provides a simple mechanism to model the dynamics of prices. Along 
similar lines it has been worked out, for instance, an option-pricing model which is more realistic than the celebrated 
Black-Scholes onej^l (recovered as the q = 1 particular case). 

The ^-exponential character of the solutions of Eq. (|14|l is not exclusive of the steady state but it also emerges 
along the time evolution 18, 23]. In the presence of multiplicative noise, the system variables directly couple to noise. 
Therefore, behaviors are observed that can not occur in the presence of additive noise alone. On the other hand 
the additive noise plays a fundamental role allowing the existence of a reasonable and normalizable steady state by 
avoiding collapse of the distribution at the origin. The particular interplay between additive and multiplicative noises, 
as well as that between deterministic and stochastic drifts, can lead to the appearance of g-exponential forms. 

The g-exponential distributions include the Boltzmann-Gibbs one as a special case (q = 1). While the latter 
corresponds to the standard thermal equilibrium, the q ^ 1 is expected to be related to quasi-stationary or long-living 
out-of-equilibrium regimes. 



B. Stochastic processes with varying intensive parameters 



Intricate dynamical behavior is a common feature of many non-equilibrium systems, which can be also charac- 
terized by power-law probability density functions. To this class belong systems whose dynamical behavior shows 
spatio/tcmporal fluctuations of some intensive quantity. This quantity may be the inverse temperature, like in the 
case of interactions of hadrons from cosmic rays in emulsion chambers; the energy dissipation in turbulent fluids, or 
simply the width of some white noise, as assumed in many financial models, such as in the famous Heston model 24]. 
The connection between this sort of dynamics and nonextensive entropy was first made by G. Wilk and Z. Wlodarczyk 
[25l | and later extended by C. Beck and E.G.D. Cohen 26], who called it superstatistics. In this "statistics of statis- 
tics" , Beck and Cohen aimed to treat non-equilibrium systems from the point of view of long-living stationary states 
characterized by a temporally or spatially fluctuating intensive parameter. Such condition can be mathematically 
expressed by 

/>oo 

B[E{z)]= /C9)e-^Wd/9 (16) 
Jo 

where B[E{z)\ is a kind of effective Boltzmann factor, E{z) a function of some relevant variable z, and f(f3) the 
probability density function of the inverse temperature /?. Superstatistics is intimately connected with nonextensive 
statistical mechanics. More precisely, it is possible to derive a generalized Boltzmann factor which is exactly B[E], 
when f((3) is the Gamma distribution, i.e., 

J bT[c]\c) 

where the g-exponential functional form of the effective Boltzmann factor turns clearly visible its asymptotic power- 
law behavior. It is noteworthy that the above effective Boltzmann factor is also a good approximation for other /(/3) 
probability density functions j2(J . 



IV. APPLICATIONS TO FINANCIAL OBSERVABLES 



A. ARCH(l) and GARCH(1, 1) processes from a nonextensive perspective 



The fluctuating character of volatility in financial markets has been considered, since a few decades ago, as major 
responsible for price change dynamics |27| . In fact, the intermittent character of return time series is usually associated 



with localized bursts in volatility and thus called volatility clustering ]28^ . The temporal evolution of the second-order 
moment, known as heteroskedasticity \2fjJ . has proven to be of extreme importance in order to define better performing 
option-price models [US 

The first proposal aiming to modelize and analyze economical time series with time- varying volatility was made by 
R.F. Engle [29j, who defined the autoregressive conditional heteroskedasticity (ARCH) process. In his seminal article, 
Engle stated that a heteroskedastic observable z (e.g., the return) would be defined as 

Z t = a t co t , (17) 

where uj t represents an independent and identically distributed stochastic process with null mean and unitary variance 
((ojt) = 0, (yit) = 1) associated to a probability density function P n (uj), and cr t the volatility. In the same work, 
Engle also suggested a simple dynamics for volatilities, a linear dependence of of on the n previous values of \zf\ , 

n 

u\ = a + J2 b i z t-i> (a,bi>0), (18) 

i=l 

afterwards named as ARCH(n) linear process f/?l /. The ARCH(n) process is uncorrelated and for n = 1 it presents 
a volatility autocorrelation function with a characteristic time of order | log &i | 33] . 

In order to give a more flexible structure to the functional form of erf, T. Bollerslev generalized Eq. I|18|) defining 
the GARCH(n,m) process^ 

n s 

a t = a + b * z t-i + c * CT * 2 -» ( c * > °) ' ( 19 ) 

i=l i=l 

which reduces to ARCH(n) process, when c 2 ; = 0, Vj. 

For the GARCH(\, 1) (pi = b and c\ = c), we can straightforwardly determine the k th moment of the stationary 
probability density function P (z), particularly its second moment 



and the fourth moment, which equals the kurtosis (k x 



«> = *. = k (i + » \_/:^_ b ^ ) . (■> + * + *k < i) . 

for processes with unitary variance, i.e., d 2 = 1. 

Continuous approaches are becoming more often used (mainly in the treatment of high-frequency data: see, e.g.. l3ll 
Isfij). Moreover, ARCH-\ike processes fail in reproducing the volatility autocorrelation function power law behavior |3fij. 
But, despite these facts, the ARCH family of processes (particularly ARCH (1) and GARCH (1, 1)) is still considered 
a cornerstone in econometrics due to its simplicity and satisfactory efficiency in financial time series mimicry. 

Having a glance at Eq. Ijl7|l . we can verify that the distribution P(z) of the stochastic variable z has, at each time step 
t, the same functional form of the noise distribution, P(u>), but with a variance at- This property allows one to look 
at process {z} as a process similar to those occuring in some non-equilibrium systems with a longstanding statio nary 
state. Specifically, this principle has allowed to establish, firstly for ARCH(1) [37J and then for GARCH(1, 1) |38| . 
a connection between b and c, P(z) and, P n (ui), the latter assumed to be of the following q„-Gaussian form 




FIG. 2: Typical runs of ARCH (1) (Upper plots) and G ARCH (1,1) (Lower plots) with a 2 = 1. For each plot the 
symbols represent the relative frequency, F(z), and the line the corresponding probability function, J^gP(z) dz. In 
(a) q = 1.2424 (6 = 0.4; q n = 1); (b) q = 1.3464 (6 = 0.4; g» = 1.2); (c) q = 1.38 (b = 0.4; c = 0.4; q n = 1); (d) <? = 1.35 
(& = 0.3;c = 0.45;<?„ = 1.2). 



By making the ansatz P (z) ~ p (z), where p (z) is the g-Gaussian probability density function which maximizes S q 
(Eq. (J3J), and by imposing the matching of second (a 2 — 1) and fourth order moments, it is possible to establish, for 
GARCH (1,1), a relation containing the dynamical parameters b and c and entropic indexes q and q n : 



6(5 - 3<z„)(2 -q) = y/{q - q n ) - [(5 - 3 9n )(2 - g) - c 2 (5 - 3g)(2 - q n )\ - c(q - q n ). (21) 

For c = 0, Eq. i|21|) reduces to the one corresponding to ARCH (1), 

= g n + 26 2 (5 - 3g n ) 
9 l + 6 2 (5-3 9n ) 1 ' 

and, for 6 = c = 0, one has q = q n . The validity of Eqs. l|2T|) and is depicted in Fig. El The discrepancy between 
p (z) and P (z) can be evaluated by computing the sixth-order moment percentual difference, which is never greater 
than 3% [2llH3. 

Since tot = Zt/cr t and (uit o t ) = 0, for q n = 1 we can write 

p (z\o * - ■ 



V2 7TCT 2 

as the conditional probability density function of z given a 2 . Considering that, 

poo 

p(z)= p(z\a 2 )P a (a 2 )d(a 2 ), 
Jo 

and P (z) ~ p (z), we obtain the stationary probability density function for squared volatility [ 
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FIG. 3: The symbols in black represent the inverse cumulative frequency, C (a 2 ) , numerically obtained for a Gaussian noise with 
b = c — 0.4 and the gray line the respective inverse cumulative distribution, f°% P a (cr 2 ) d (cr 2 ) with (k, A, <r 2 ) = (0.444, 2.125, 1) 
for Pa (a 2 ) . 

where 

A- 1 1 k- 1_g 
q-1 2' o- 2 (3q- 5) 

As one can observe in Fig. the ansatz gives a quite satisfactory description for a 2 probability density function, 
suggesting a connection between the ARCH class of processes and nonextensive entropy. These explicit formulas 
can be helpful in applications related, among others, to option prices, where volatility forecasting plays a particularly 
important role. 

Albeit uncorrelated, stochastic variables {z t } are not independent. Applying the g-generalized Kullback-Leibler 
relative entropy |39, 40j to stationary joint probability density function p\ (z t ,z t -i) andp2 (z t ,z t -i) = p (zt) p (zt-i) = 
[p (z)] it is possible to quantify the degree of dependence between successive returns, through an optimal entropic 
index, q op . In Ref. [38j. it was verified the existence of a direct relation between dependence, q op , the non-Gaussianity, 
q, and the nature of the noise, q n . An interesting property emerged, namely that, whatever be the pair (b, c) that 
results in a certain q for the stationary probability density function, one obtains the same q op and, consequently, the 
time series will present the same degree of dependence 38j . See Fig. 

It was also verified (see Fig- EJ) for ARCH (1) that the degree of dependence varies visibly with b and with the lag r 
between returns. This dependence would be related to a short-memory in volatility j^jj . The variance between these 
results and the empirical evidence of persistence in the real return time series dependence degree for time intervals up 
to 100 days recently found, shows that financial markets dynamics are, in fact, ruled by some form of long-memory 
processes in volatility 01 ■ 

Let us comment at this point that the connection between several entropic indexes, each one related to a different 
observable, is compatible with the dynamical scenario within which the nonextensive statistical mechanics is for- 
mulated. In fact, several entropic indexes emerge therein, coalescing onto the same value q = 1 when ergodicity is 
verified. 
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FIG. 4: Plot of q op versus q for typical (b,c,q n ) triplets. The arrow points two examples which were obtained from different 
triplets and nevertheless coincide in what concerns the resulting point (q,q op ). 




FIG. 5: In panel (a) q op versus b for ARCH (1) process. Even for large allowed b values, the decrease of dependence degree 
(i.e., increase of q° p ) is visible when the time lag increases, which is not compatible with the empirical evidence of dependence 
degree around q° p — 1.275 verified in Dow Jones return time series for time lags up to 100 days, panel (b). Comparing the two 
figures (a) and (b), it appears that financial markets present long-memory in the volatility similar to a ARCH (n) process with 
n > 1. 



B. Mesoscopic models for traded volume in stock markets 



Another important observable in financial markets is the traded volume V (the number of shares of a stock traded 
in a period of time T). In Ref. 42] it was proposed an ansatz for the probability density function of high-frequency 
traded volume, which presents two power-law regimes (small and large values of V), 

P (v) = ^T^(-l)> (23) 

where v represents the traded volume expressed in its mean value unit (V), i.e., v = V/(V), a and 9 are positive 
parameters and Z = ^{v/9)~ a exp {—%) dv. 

The probability density function (|23|) was recently obtained from a mesoscopic dynamical scenario based in the 



ulity 

Feller process [44| (using Ito definition) 



<lr - [v- J dt+^2vl dW t , (24) 



where instead of being constant in time, f3 varies on a time scale much larger than the time scale of order 7 -1 needed 
to reach stationarity. The deterministic term of Eq. 124J1 represents a restoring market mechanism which tries to keep 
the traded volume in some typical value = ^rp- and the second term reflects stochastic memory and, basically, the 
effect of large traded volumes. In fact, large values of v will provoke large amplitude of the stochastic term, leading 
to an increase or decrease of the traded volume (stirred or hush stock) depending on the sign of Wt . The fluctuation 
of (3, alike to fluctuations in the mean value of v, can be related with changes in the activity volume due to agents 
herding mechanism caused by price movements or news. 

Solving the corresponding Fokker-Planck equation for Eq. I|24|) we got the conditional probability of v given /3, 



P(v\0) 



P 



{f3v) a exp (-/?«) (a > -1, /3 >0). 



(25) 



Since 



P(v) 



P{v,(3) dp 



p(v\/3) P(/3) d/3, 



(26) 



considering that (3 follows a Gamma distribution, 



P(P) 



1 



XT [5] \\ 



s-i 



exp ( -- 



Eq. |2E|) yields, 



. . I fv\ a 



(S > 0, A > 0), 



(-!)■ 



(27) 



(28) 



where, q = 1 + 1/(1 + a + S) and 9 = (q — 1)A. A numerical simulation for 1 minute traded volume in NYSE stock 
index is exhibited in Fig. 
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FIG. 6: In panel (a) open symbols represents the PDF for the ten-high 1 minute traded volume stocks in NYSE exchange; solid 
symbols represent the PDF obtained for the numerical realization depicted in panel (b) and line the theoretical PDF Eq. 1281 1. 
Parameters are q = 1.17, a = 1.79, A = 1.42 and S = 3.09. 

Another possible mechanism to describe the dynamics of volumes, has been recently proposed^5| on the basis of 
mean-reverting processes with additive-multiplicative structure, namely, 

dv — —7 (v — 9) - dt + nVv dWi + adW 2 , (29) 
v 

where, a, [i are positive constants and W\, W2 independent standard Wienner processes. That is, following the 
ideas presented in Section fill Al in addition to fluctuations endogenously processed by the market, other fluctuations 
are taken into account that affect the dynamics directly. The stationary solution of the Fokker-Planck equation 
associated to the Ito-Langevin Eq. H29fl has the form H28[). Moreover, the underlying sequences present intermittent 
bursts (similar to those exhibited in Fig. |HJb) and, preliminary results indicate the presence of persistent power-law 
correlations, as observed in real data sequences. 

Eq. H29JI belongs to a larger class of processes with two-fold power-law behavior that can be also suitable for 
volumes, as well as, for other mean-reverting variables such as volatilities. 



V. FINAL REMARKS 



Additive-multiplicative processes are at the core of nonextensive statistical mechanics. In the same natural way that 
standard Brownian motion leads to Gaussianity, linear additive-multiplicative stochastic processes lead to g-Gaussian 
distributions. Special (scale-invariant) correlations, that forbid convergence to the usual Gauss or Levy limits, lead 
to a new type of statistical distributions. A remarkable feature is that the resulting power-law distributions may have 
exponents out of the Levy range, thus allowing to embrace a larger variety of empirical processes. The presence of 
two Gaussian white noises, one either enhanced or reduced by internal information, and another purely exogenous, 
represents quite realistic features present in a variety of systems, thus justifying the ubiquity of the probability 
distributions associated to such kind of processes. In particular, as we have shown, they allow to model the dynamics 
of prizes, volatilities, stock-volumes and other relevant financial observables. 

Another expression of the mechanism underlying the nonextensive theory is connected to the existence of a fluc- 
tuating intensive parameter (or "inverse temperature" ) following the ideas that foundate the Beck and Cohen super- 
statistics [2^. We have shown that these principles allow an alternative description of the dynamics of stock- volumes. 
Furthermore, such kind of mechanism allows an interesting perspective for treating the family of G/ARCH processes. 
The fact that the resulting probability density functions can be described in terms of g-Gaussian distributions, pro- 
vides a tractable way of dealing with empirical distributions that match the G/ARCH types. Some of the discussions 
presented in this review have been done at a mesoscopic scale. The determination, from more microscopic models, of 
the parameters used at the mesoscopic scale is certainly welcome. 
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